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Abstract 

The correction to the muon anomalous magnetic moment from the pion-pole contribution to the 
hadronic light-by- light scattering is considered using a description of the tt° — 7* — 7* transition 
form factor based on the large- Ac and short- distance properties of QCD. The resulting two- 
loop integrals are treated by first performing the angular integration analytically, using the 
method of Gegenbauer polynomials, followed by a numerical evaluation of the remaining two- 
dimensional integration over the moduli of the Euclidean loop momenta. The value obtained, 
a L ^° = +5.8 (1.0) x 10~ 10 , disa grees with other recent calculations. In the case of the 
vector meson dominance form factor, the result obtained by following the same procedure reads 
a^ byL;7r ° I vmd = +5.6 x 10~ 10 , and differs only by its overall sign from the value obtained by 
previous authors. The inclusion of the rj and 7/ poles gives a total value aJf yL;PS = +8.3 (1.2) x 
10 -10 for the three pseudoscalar states. This result substantially reduces the difference between 
the experimental value of a M and its theoretical counterpart in the standard model. 
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1 Introduction 



A high-precision measurement of the muon anomalous magnetic moment is presently being 
performed at the Brookhaven National Laboratory by the Muon (g — 2) Collaboration The 
value released recently [||, 



a M+ = 11659 202 (14) (6) x 1(T 10 , (1.1) 

already improves the precision by a factor of 8 as compared to the previous determination 
at CERN 0. The final aim of the experiment is to reach a precision of ±4 x 10~ 10 , after 
completion of the analysis including the full set of data collected during recent years. 

If only contributions of (multiflavored) QED are taken into account, the theoretical accuracy 
is still much better than the present experimental one [|J [|. However, hadronic effects, which 
certainly cannot be ignored in the case of the muon g — 2 factor, are unfortunately more 
difficult to estimate, and the corresponding uncertainties substantially increase, and actually 
completely dominate, the error of the total theoretical value. Nevertheless, the value (|1 . 1|) 
disagrees with various theoretical calculations in a way that could become very significant as 
the BNL E821 experiment further decreases its statistical error. Seen from this perspective, it 
is not surprising that the announcement of the value ( |1.1| ) has triggered a wealth of theoretical 
activity (a complete bibliography of the recent work in this direction would represent a tedious 
task ||). It is, however, our feeling that, rather than calling for hasty conclusions concerning 
possible evidence for physics beyond the standard model, the present situation first requires 
that the hadronic contributions be scrutinized anew and with greater care. 

The hadronic contributions to the anomalous magnetic moment of the muon considered 
so far naturally fall into three categories: vacuum polarization, light-by-light scattering, and 
higher-order electroweak contributions, as represented in Fig. |T[ Hadronic effects in two-loop 
electroweak corrections are estimated in Refs. |7], || ||. They are small, of the order of the 
expected experimental error, and the associated theoretical uncertainties, of the order of 10%, 
can be brought under safe control 0. Hadronic vacuum polarization effects benefit from the 
fact that they can be related to the total cross section for e + e~ — > hadrons [TO]. However, 



existing data come from different sources, and do not always meet the required accuracy, so 
that they are supplemented with theoretical input. There exist therefore different estimates 
for the hadronic vacuum polarization, which either strengthen the difference between theory 
and the experimental value (|1.1|) , or make it appear only marginal. Our purpose here is not 
to intervene in this ongoing debate [11, 12, 13, 14, [15|, 16, T7|, but rather to have a closer look 



at the remaining item on the list, namely, the hadronic effects in the so-called light-by-light 
scattering contribution. There exist several estimates of the latter [IS, |19|, |2(], 21, |22j. The 



latest to date, obtained by two different groups (see also j23|), give values that are consistent 
within the quoted errors, 

a LbyL;had = _ 7. 0. ( 1 . 5 ) X 1 () ~ 10 [20,13, ( L2 ) 
a LbyL;had = _C; . 2 (3 . 2 ) X 1 " 10 |2J. (1.3) 

On the other hand, both groups have changed their value of a^ yU had , and even its sign, at 
some stage, and for reasons that are not very easy to follow - partly because both analyses 
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Figure 1: The three topologies involving hadronic contributions to the anomalous magnetic 
moment of the muon: from left to right, vacuum polarization insertion in the vertex, light-by- 
light scattering, and two-loop electroweak contributions. The cross indicates an insertion of the 
electromagnetic current, the shaded areas correspond to hadronic subgraphs, while the lines 
exchanged in the rightmost graph correspond to a photon and a neutral gauge boson. 



rapidly leave the realm of analytical work and resort to numerical techniques. This aspect of 
the present theoretical situation provided the motivation to improve from the analytical side 
the study of the hadronic light-by- light contributions. 

A complete discussion of the hadronic light-by-light contributions involves the full rank-four 
hadronic vacuum polarization tensor Hn V \p(<li, 92 > Q3), which is a rather involved object. On the 
other hand, the values ( |1.2| ) and ( |1.3| ) are both dominated by a well- identified component, 

~LbyL; had -^LbyLjTr^ _j_ 

u p _ M <"''■> 

aJf yL; "° = -5.7 (0.3) x 1(T 10 [21] , (1.4) 
a™ = -5.9 (0.9) x 10- 10 @ , (1.5) 

arising from the one-particle reducible pion-exchange pieces of H^XpiQi, Q2, 93)- In the present 
paper, we focus on the latter, leaving for future work a complete discussion of H^uXpiQi, q 2 , q^) 
within a theoretical framework close to the one adopted here (and discussed below) for the 
pion-pole contribution. 

The main ingredient in the determination of the pion-exchange graphs depicted in Fig. |2| is 
the double off-shell pion-photon-photon transition form factor F^y*-/* (Qi, No sufficiently 
complete experimental information on this quantity is available at the time being, so that even 
the one-pion-exchange component of Hp,i>x P {qi, q 2 , (fe) is beyond reach from an experimental 
point of view. Resort to theory is therefore unavoidable. Now, from the QCD perspective, the 
theoretical knowledge of T^^y* is sparse. There exists a well-defined limit of QCD where the 
situation improves somewhat, namely, the limit of infinite number of colors Nc @, j25|. As Nc 



becomes large, the spectrum of QCD reduces, in each channel with given quantum numbers, 
to an infinite tower of zero-width resonances. As a consequence, the analytic structure of the 
form factor T^o^^i^q^q^) becomes simpler: it consists of a succession of simple poles, due to 
the contributions of zero- width J PC = 1 states (e.g., the p meson and its radial excitations) 
in each channel. Dealing with an infinite number of resonances remains cumbersome, and to 
a large extent illusory, since the characteristics of these states (masses, couplings, etc.) are 
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Figure 2: The pion-pole contributions to light-by-light scattering. The shaded blobs represent 
the form factor J r 7r o 7 * 7 *. The first and second graphs give rise to identical contributions, in- 
volving the function T 1 (g 1 ,g 2 ;p) in Eq. ( |3.4| ), whereas the third graph gives the contribution 
involving T 2 (q 1 ,q 2 ;p). 



in general not known. It has, however, been shown in several instances (see, e.g., the review 
F2B | and references therein) that keeping, in each channel, only a finite number of resonances, 
supplemented with information on the QCD short-distance properties coming from the operator 
product expansion |27|, |28j, already gives a good description of quantities like form factors or 
correlation functions in the Euclidean region, especially when they occur in weighted integrals 
over the whole range of momenta. In particular, the form factor T^o^*^* has recently been 
studied p!| from the point of view of this lowest meson dominance (LMD) or minimal hadronic 



Ansatz (MHA) approximation to large- Nc QCD. Since the analyses carried out in Refs. [^T], ^| 
rely on models of jF 7r o 7 » 7 », vector meson dominance (VMD) or extended Nambu-Jona-Lasinio 
(ENJL) (see [22 and references therein), that do not reproduce the correct QCD short-distance 
properties, a second motivation for the present study was to compare the results ( |1.4| ) and ( |1.5| ) 
with those derived from a representation of J r 7r o 7 . 7 . that complies with these QCD constraints. 
Finally, let us mention for completeness that the pion-pole contribution we are interested in 
corresponds to the lowest-mass part of Tl^uXpiqi, ?2, #3) that is leading in the large- Nc limit f3~0fl , 
which might provide an explanation as to why it happens to constitute the dominant fraction 
of the light-by-light scattering correction to a^. 

The remaining material of this paper is organized as follows. Section recalls a few defini- 
tions that are relevant for the light-by-light contribution to a^, and also serves the purpose of 
introducing our notation. The expression for the two-loop integral which gives the pion-pole 
contribution a^^ to the anomalous magnetic moment of the muon in terms of jF 7r o 7 » 7 » is de- 
rived in Sec. 0. In Sec. (|we discuss a generic class of form factors to which those inspired from 
large- Nc QCD and considered here belong, but that also covers other cases, like the constant 
form factor, given by the Wess-Zumino-Witten term []31| , or the vector meson dominance form 
factor. The method of Gegenbauer polynomials is presented in Sec. [5] and used in order to 
perform the angular integrations. This then leads to a two-dimensional integral representation 
for a^^ in terms of .F 7r o 7 * 7 « and of several weight functions (Sec. |]). Numerical results for 
a LbyL;7r° are presented in Sec. [7], where the contributions from the r\ and rf poles are also briefly 
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discussed. A discussion and conclusions make up the content of Sec. ||. For the reader's conve- 
nience, the properties of the form factor jF 7r o 7 * 7 , that form the basis of the representations we 
use here have been gathered in an Appendix. 



2 Definitions 

The interaction of the photon field with the standard model fermionic degrees of freedom is 
described by the Lagrangian (e denotes the electric charge of the electron) 



£(^7 P ^)W " E e <^7 p <?)( 



x 



(2.1) 



The first sum runs over the charged lepton flavors i = e~ fi~ r~, while the second sum runs over 
the quark flavors, with e q standing for the corresponding electric charges expressed in units of 
|e|. Our interest hereafter is restricted to the contributions of the light quarks q = u,d, s. The 
muon gyromagnetic ratio is obtained from the muon proper vertex function T p (p',p), defined 
as (p 2 = p 12 = m 2 , m denotes the muon mass) 



m(p')T p (p',p)u(p) = (rt')IHe)£(^7p^)(0) + ( ie ) £ e q (q lp q)(0) ^'(p)) , (2.2) 



at vanishing momentum transfer. 

The contribution to F p (p',p) of relevance here is the matrix element, at lowest nonvanishing 
order in the fine structure constant a = e 2 /(4n), of the light quark electromagnetic current 



2 1 1 

3p(x) = ^lpu){x) - -{d^ p d){x) - ~(s7 p s)(x) (2.3) 



between \i states, 



-ie)u(p')f p (p',p)u(p) = (// (p')\(ie)j p (0)\fi (p)} 

d 4 q 1 f d 4 q 2 (-i) 3 



x 



(2tt)*J (2tt)4 qiq 2 2 (qi + q2-k) 2 
i i 



(p' — qi) 2 — m? (p' — qi — q^) 2 — m 2 
x(-te) 3 u{p')Y{p"- di+mWip"- 4i- ^ 2 + m) 7 A u(p) 
x (ze) 4 n^ Ap (gi, q 2 , k - q x - q 2 ) , (2.4) 



with k p = {p' — p) p and 



IW(?i,<72,<73) = / d A x x J d 4 x 2 J d*x 3 e i( -^+^ + ^ 

x(n\ Tij^xxKvMhMjpio)} | n > (2.5) 

the fourth-rank light quark hadronic vacuum polarization tensor, | Q ) denoting the QCD vac- 
uum. 
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Since the flavor diagonal current j^x) is conserved, the tensor n M „> p (gi, q 2 , qs) satisfies the 
Ward identities 

{q.x\ Q3 ; (?i + Q2 + g 3 ) p }n^ Ap (gi, 92, *$) = . 

This entails that Il^ Ap (gi, q 2 , k-q 1 -q 2 ) = -fc ff (0/0fc p )n /u/Aff (gi, g 2 , k - q 1 - q 2 ) andQ 



(2.6) 



u(p')f P (y,p)u(p) = u (p0 

as well as r p (p',p) = k a T pa .(p',p) with 
u(p')fW(p',p)u(p) = 



7p F 1 (A; 2 ) + -a pT rF 2 (fc 2 ; 
zm 



u(p) 



(2.7) 



x 



d i q 1 f d A q 2 



1 



(2^7 (2vr)Ni 2 g 2 2 (?i + 92-fc) 2 
1 1 



(P' — <?l) 2 — m2 (P' — <?1 ~ ?2) 2 — m2 

xu{p')^(f- ^i+m) 7 I '(y / - ^1- ^ 2 + m) 7 A u(p) 
.9 

x ^uXaiQl, <?2, k-qi-q 2 ). 



{21 



Following Ref. |32| and using the property k p k a u(p')T pa (p' , p)u(p) = 0, one deduces that 
-Fi(O) = and that the hadronic light-by-light contribution to the muon anomalous magnetic 
moment is equal to 



(0) = -^-tr{(^ + m)[ 7 ", 7 CT ](^+m)f p(J (p,p)} . 



48m 

Equivalently, one may project out F 2 {k 2 ) from T p (p',p), 

F 2 {k 2 ) = tr {(p' + m) A p 2 \p',p)(p" +m)T p (p',p)} 



with the help of the projector |33| 
A^(p',p) 



m 



k 2 (4m 2 - k 2 ) 



k 2 + 2m 2 , ' 
7 " + m(k 2 - Am 2 ) {P +P)p . 



(2.9) 



(2.10) 



(2.11) 



One then uses T p {p',p) = k a T p ° (p' ',p) and employs the identity (for p 2 = p' 2 = m 2 ) 



+ m) 7 P (p" + m) = {p 1 + m) 



(p" + m) (2.12) 



to simplify A p 2 \p' ,p). Averaging over the directions of the four- vector k^ with the constraints 
k ■ p = —k 2 /2, k ■ p' = k 2 /2 and keeping the leading terms in k in the resulting expression then 
gives back Eq. (p|) E§ B3. 



1 We use the following conventions for Dirac's 7 matrices: {7^, 7^} = 2?7 M „, with 7y M „ the flat Minkowski space 
metric of signature (+ — — — ), a^ v = 7^], 75 = i7°7 1 7 2 7 3 , whereas the totally antisymmetric tensor 

s^ppa is chosen such that £0123 = +1- 
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3 Pion-pole contribution 

From now on, we concentrate on the contributions to H^\ p {qi, g 2 , Q3) arising from single neutral 
pion exchanges (see Fig. 0), which read 



2^ 



n^ Ap (9i> 92, g 3 ) = 1 — 21:1 ( gi + g 2 ] 2 _ M 2 <& e ^ ^3 (gi + 52) 

. ^vo 7 «7*(gi' (gi + g2 + gs) 2 ) ^o r7 *(g|, gf) „, ^ j t 

( g2 + g 3 )2 - M 2 MPQ/3 ?1 ^ 2 *^ ?2 % 

, pTfirAlh ?j) ^rVrjgl (gj + g2 + g 3 ) 2 ) « /3 ^ r 

+ * ^ + g 3 )2 _ ^2 £/xA«/9 gl g3 £ *>p<rr <? 2 Igl + g3j ■ 

(3.1) 

The form factor JT 7r o 7 » 7 »(q' 2 , gf) is defined as (using the same convention as in Ref. p9| ) 

i I d 4 xe^(n\T{j,(x)UO)}k°(p)) = WgVJvv 7 *(g 2 .(P-g) 2 ) ( 3 - 2 ) 



with jF 7r o 7 . 7 . (g 2 , gf) = jF 7r o 7 * 7 » (g|, g 2 ). For the computation of the corresponding value of 

a LbyL;.« = F 2 (0 ) | pion polo , we need 



An 



n^2(g 1 ,g 2 ^-g 1 -g 2 ) 



J 7 n .o 7 * 7 * (g 2 , 

g|) ^o 7 * 7 »((g! + g 2 ) 2 ,0) 
(gx + g 2 ) 2 -M 2 
^V°7*7* (gi , 0) ^ r 7r o 7 . 7 . (gf , (gi + g 2 ) 2 ^ 



(gl + 92) 1 



~r 4 2 S/f2 ^parpHl ^v\ol^H\H2 

gl ~~ iW 7T 

. -7vy 7 * (g 2 ; (gi + g2) 2 ) Tj^fj* (gf , 0) Q p T 

~i~ 1 2 A/f2 £ p\a/3 Qi Q2 e vapr Q2 

g 2 ~ IVI TT 

+ 0(k). (3.3) 

Inserting this last expression into Eq. ( |2.8| ) and computing the corresponding Dirac traces (we 
used REDUCE ||) in Eq. (U), we obtain 



= -e 6 



cfgi r drq 2 



^ (0«r\4 / (O^M n 2 n 2 



with 



(2tt) 4 J (2tt) 4 g 2 g 2 2 (gi + g 2 ) 2 [(p + gi) 2 - m 2 ] [(p - g 2 ) 2 - m 2 

(gi 2 , (gi + g2) 2 ) J>o 7 . 7 . (gf , 0) 
x 2 — ^t 2 J Hgi,g2;pj 

g2 - ^ 

^07*7* (<? 2 , gf) ^vo 7 . 7 * ((gi + g 2 ) 2 , 0) 



16 16 
^i(gi, g2!P) = y (p • gi) (p • g2) (gi • to) - -j (p • ?2) 2 g 2 



(3.4) 



- 1 (p • 9i) (gi ■ 92) g 2 + 8 (p • 92) g 2 g 2 , - y (p ■ 92) (gi • g2) 2 
16 16 

+ — m 2 qlql - — m 2 (q 1 ■ q 2 ) 2 , (3.5) 

16 16 

72(gi, q.2\V) = y(p- ?i) (p • ?a) (?i -92) - y (p • 9i) 2 92 

+ g (P ■ ?i) (9i ■ 92) 92 + g (P ■ 9i) 9i 92 

+ g ™ 9l 9 2 - g 111 (9l • 92) • (3.6) 

In deriving Eq. ( |3.4j ), we have used the fact that, upon a trivial change of variables in the two- 
loop integral (|2.8|) , the two first terms of Eq. ( |3.3|) lead to identical contributions. Furthermore, 
in writing T 2 (gi, q 2 ; p) we have taken into account the invariance of the remaining factors of the 
corresponding integrand under the exchange q\ <-> — q 2 - 

Before discussing the pion-photon-photon transition form factor, let us briefly mention a 
few features of Eq. (|3.4j) that are independent of the precise form of T^*^* . One observes that 
there are five independent variables p ■ qi,p- 92, 9i> gf, an d 91-92, whereby p ■ qt only occur in the 
fermion propagators, not in the form factors. It might therefore be possible to perform these 
two integrations over p-q\ and p-q 2 for general form factors. We would then be left with a three- 
dimensional integral, since the form factors depend on q\ ■ q 2 . Note that in both contributions 
in Eq. (|3.4J) there is one factor jF 7r o 7 « 7 « (g 2 , 0). In principle, this function could be extracted, 
over a finite energy range, from the data collected, e.g., by the CLEO collaboration |37|]. This 
would reduce the model dependence. 



4 The pion-photon-photon transition form factor 

In the present work, we shall consider a representation of the form factor jF„-o 7 * 7 *, based on 
the large- -/V^ approximation to QCD, that takes into account constraints from chiral symmetry 
at low energies, and from the operator product expansion at short distances. Our numerical 
estimates for g — 2 below will be mainly based on expressions for the form factor jF 7r o 7 * 7 * that 
involve either one vector resonance (lowest meson dominance, LMD) or two vector resonances 
(LMD+V). These Ansdtze have been thoroughly discussed in Ref. |29[, and a short summary 
can be found in the Appendix. Furthermore, as a reference point we shall also consider the 
simplest model for the form factor that follows from the Wess-Zumino-Witten (WZW) term 
that describes the Adler-Bell-Jackiw anomaly in chiral perturbation theory. Since in this 



case the form factor is constant, one needs an ultraviolet cutoff, at least in the contribution to 
Eq. ( |3.4j ) involving T x . Therefore, this model cannot be used for a reliable estimate, but serves 
only illustrative purposes and as a check of our calculation. Finally we shall use the usual 
vector meson dominance form factor. Essentially, this corresponds to the form factor that 
has been employed in previous calculations for the pion-exchange contribution to g — 2 of the 
muon |19|, |21|] ; see also Ref. [ 3J| . This will allow us to check the numerics. The expressions 



for the form factor jF 7r o 7 * 7 * based on the ENJL model that have been used in Ref. J22] do not 



allow a straightforward analytical calculation of the loop integrals. However, compared with 
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the results obtained in Refs. p0| , 21] , the corresponding numerical estimates are rather close 
to the VMD case (within the error attributed to the model dependence). We would like to 
stress again that both the VMD and ENJL form factors fail to correctly reproduce the QCD 
short- distance constraints discussed in As noted earlier, the question of how sensitive the 
results are to the latter motivated this work. 

For the four cases mentioned above, the form factors are given by |29|, |39]] [different choices 
for the global sign in the normalization of the form factors are of no relevance, cf. Eq. ([3.41)1 



T WZW( 2 2n 



T VMD ( 2 n 2\ 

j- 7r o 7 . 7 4g 1 , q 2 ) 



T LMD l 2 2\ 

J- 7r o 7 * 7 4g 1 , q 2 ) 



with (see the Appendix) 



N, 



c 



12tt 2 F w 

N c 



Ml 



Ml 



12tt 2 F w {ql 



Ml) 
ql + q 2 2 -cy 



My) 



3 

^ 5?«§(9i + 9: 
3 



(ql - Ml){q\ - Ml) 



HQi + 4? + + h(ql + q 2 2 ) + /l 7 



(g 2 — M Vl ) (ql 



Ml 2 )(ql 



MlM 



Ml 2 ) 



cv 



4tt 2 Fl ' 



/l7 



N c M Vl M* 2 



4tt 2 



(4.1) 
(4.2) 
(4.3) 
(4.4) 



(4.5) 



According to Ref. [fiO| |, the form factor !F n o 7 * 7 * (— Q 2 , 0) with one photon on shell behaves 
like 1/Q 2 for large spacelike momenta, Q 2 = —q 2 . The experimental data are compatible with 
this behavior |3"7fl . Whereas the LMD form factor does not have such a behavior, it can be 
reproduced with the LMD+V Ansatz, provided that hi = 0. In Ref. 
form factor to the CLEO data yielded h x = -0.01 ± 0.16 GeV 2 , h 5 



29 



, a fit of the LMD+V 
6.88 ± 0.61 GeV 4 , or, 

if we fit /15 with hi = 0, h 5 = 6.93 ± 0.26 GeV 4 . On the other hand, the parameter h 2 that 
appears in Eq. (4^) is not known. If we compare, in the context of Ref. the expression 
for the 7T° — > e + e~ partial width based on the LMD+V form factor to the experimental values 

h 2 \ ^20 GeV 2 , although slightly positive 



42, 43 



we obtain a rather loose constraint on h 2 , 
values for h 2 seem to be preferred in order to reproduce the experimental rate. 

A crucial observation at this stage is that all dependences on qi ■ q 2 in the numerators in 



lmd+v,lmd^2^ (q 1 -\-q 2 ) 2 ) can be canceled by the ones in the resonance propagators. Therefore, 
we can write 



7r°7*'7 



J- n-o 7 * 7 * (q-^ , q 2 ) = 



fill) - E 



qi - M, 



-9m V4 (ql) 



v, 



(4.6) 



This representation, which seems to hold quite generally in the large-A^ limit of QCD (see 
the Appendix), obviously extends to the VMD and WZW form factors. The corresponding 
functions f(q 2 ) and gu v . (q 2 ) can easily be worked out from the explicit expressions fl4.1|)-(fL4|). 
They are displayed in Table |l| For the VMD and LMD form factors, the sum in Eq. fl4.6|) 
reduces to a single term, and we call the corresponding function gM v (q 2 )- In the case of the 
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Table 1: The functions f(q 2 ), g Mv (q 2 ), and g(q 2 ;M 2 ) of Eqs. ([Oj|) and Q) for the different 
form factors in Eqs. ( |4.1|) -( f4~4|) . 





f(a 2 ) 

j \q ) 


n** ( n 2 \ 


wzw 


N c 
4tt 2 F 2 





VMD 





N c M v 
Atx z p 2 q l — My 


LMD 


1 


q 2 + M 2 - c v 


q 2 -M 2 


q 2 -M 2 


LMD + V 


q 2 + h 


g(q 2 ; M 2 ) = q 2 M 2 (q 2 + M 2 ) + h(q 2 + M 2 ) 2 
+ h 2 q 2 M 2 + h 5 (q 2 + M 2 ) + /i 7 


(q 2 - M^W - M 2 2 ) 



LMD+V form factor, we have written 



9m Vi {q 



1 g(q 2 ;M 2 



M 2 2 -M Vi (q 2 - M Vi )(q 2 - M 2 2 ] 



( 2, = I g(q 2 ;M 2 2 j 

9Mv > [q) M 2 -M 2 {q 2 -M 2 ){q 2 -M 2 - 



L Vi lvl v 2 vy lv± Vi)v± lvl v 2 J 

with the function g(q 2 ; M 2 ) given in the table. The expression ( |4.6| ) also holds for q\ = 0, which 
implies that /(0) = in order to reproduce the asymptotic l/Q 2 behavior of J 7 ^*^* (— Q 2 , 0). 
As noted above, this fails to be the case for the LMD Ansatz f LMD (0) ^ 0, and holds for the 
LMD+V Ansatz provided h\ = 0. Note that the limit Q\ — > oo with Q\ fixed can actually not 
be studied with the operator product expansion (see the discussion in Ref. [F2U{). 



Thus, using partial fractions, we can write the product of the form factors in the two 
contributions to the right-hand side of Eq. ( |3.4|) as follows: 



•7v>7*7* (qf, (qi + q2) 2 ) -Tyy^* (gf, 0) F n ^ r 7r o 7 * 7 * (gf, 0) 
g 2 2 -M 2 3 gf-M 2 



(9i + qi) 2 ~ M 2 J [ - 



and 



J r 7T Qj*~ f *(q 1 ,q 2 ) ^ r 7rQ7*7*((gi + ga) , 0) _^-k / 2 2^ 

(<?i + <? 2 ) 2 - M 2 " 3 /lVf(9l,?2 ' 
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X 



( gi+?2 )2_ M 2 



no) + y 9Mv ^ ] 

h. {Ml- Ml) 



+ y g J^M | (49 ) 

Note that in Eq. (}4.9|) all terms contain a propagator + q^) 2 — M 2 ], M = M n or My, in 

contrast to the prefactor of 7\ shown in Eq. ( |4.8| ). 

Since the form factors discussed above go to a constant for small momenta, the integrals 
in Eq. (|3.4j) are infrared finite. The behavior at large momenta also ensures the ultraviolet 
convergence of these integrals, except, of course, for the constant form factor T^Efl*- However, 
the integral involving T2(qi, q2',p) is finite in the latter case also (this has been observed before 



21 1, and we have confirmed it numerically as well). 



5 Angular integrations 

As we shall show in this section, for form factors ^ r 7r o 7 * 7 . (qf, q%) that have the general form 
given by Eq. ( [4.6|) , it is possible to perform all angular integrations in the two-loop integral 
of Eq. (|3.4|) using the technique of Gegenbauer polynomials (hyperspherical approach); see 
Refs. [[0], [45], |35| . In order to do so, we perform a Wick rotation of the momenta gf, q% , and p^, 
denoting by capital letters the rotated Euclidean momenta, with Qf = —qf, P 2 = — m 2 , etc. 



Let us mention that in the language of Refs. [45, 35fl the loop integrals we have to deal with are 



planar. We also note that these integrals are of the two-loop self-energy type. Other techniques 
that lead to either one- dimensional dispersive or two-dimensional integral representations [[1(J 
could a priori also be relevant in the present context, although probably not without further 
specifying the functions f(q 2 ) and gM v .(q 2 ) m Eq- 



5.1 Properties of Gegenbauer polynomials 

Let us briefly summarize some basic properties of the Gegenbauer polynomials; see also 
Refs. |48|. We write the measure of the four-dimensional sphere as follows (Euclidean 
momenta): 

d 4 K = K 3 dK dVL(k) , / dtt(K) = 2tt 2 . (5.1) 
The generating function of the Gegenbauer polynomials^ C n (x) is given by 

1 oo 

J2z n C n (x), -1 < x < 1, \z\ < 1. (5.2) 



z 1 -2xz + l n=0 

From Eq. ( |5.2| ) we immediately obtain the following property under parity transformations 
C n (—x) = (— l) n C n (x). Furthermore we get C n {l) — n + 1. The Gegenbauer polynomials obey 

2 We shall need only the special case C„(x) = C„\x). 
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the orthogonality conditions 



/ dQ(K) C n {Q x ■ K) C m (K ■ Q 2 ) = 27T 2 ^-C n (Q 1 -Q 2 ), 
J n + 1 

Jdn(k)C n (Q-k)C m (K-Q) = 2vr 2 5 nm , (5.3) 

where, for instance, Q\ ■ K is the cosine of the angle between the four-dimensional vectors Q\ 
and K. Some low-order cases of the polynomials are given by C (x) = 1, C\(x) = 2x, C 2 (x) = 
Ax 2 - 1 and therefore x = C l (x)/2, x 2 = [C 2 (x) + C (x)]/4. 

From the generating function, we obtain the following representation of the propagators in 
Euclidean space: 

i yM oo 

hL J2( Z K L ) n C n (K-L), (5.4) 



(K-L) 2 + M 2 \K\\L 



n=0 



-i yM oo 

^ L J2(-ZM L ) n C n (k.L), (5.5) 



(K + L) 2 + M 2 \K\\L\ 



n=0 



M K 2 + L 2 + M 2 - y(K 2 + L 2 + M 2 ) 2 - 4K 2 L 2 

Zkl = 2\k¥\ ' (5,6) 

Note that we have to choose the negative sign in front of the square root in Z^ L in order that 
\Zft L \ < 1. For a massless propagator these expressions simplify as follows: 

1 »(i-m±±(m'<**-L) 



(K-L) 2 V \K\J K 2 ^ \\K 

1 

T 2 



-r 2 t( l -m) c n {k-L). (5.7) 



n=0 



5.2 Basic angular integrals 

Let us introduce the following abbreviations for the propagators in the loop integral in Eq. 

D 1 = Ql, D 2 = Q 2 2 , ^ 3 = (Qi + Q 2 ) 2 , D^iP + Q.f + m 2 , 

D 5 = (P-Q 2 ) 2 + m 2 , D M = (Q 1 + Q 2 ) 2 + M 2 , (5.8) 

where M denotes either the pion mass or the mass of some vector resonance. For M — ► we 
recover the photon propagator D 3 . 

By rewriting the scalar products P-Qi, P-Q 2 , and Q\-Q 2 in terms of the propagators -D4, D 5 , 
and -D3, respectively, we can remove some of the terms in the numerator of T n j (D1D2-D3-D4-D5), 
n = 1,2, in Eq. (|3.4|) . If we multiply by the form factors, taking into account their general 
form from Eqs. ( |4.8| ) and ( |4.9| ), and use a partial fraction decomposition, we finally obtain the 
following basic angular integrals [apart from the trivial one J dQ(Qi)dQ(Q 2 ) = 47r 4 ], which 



11 



can all be performed using the method of Gegenbauer polynomials (from now on, we write Qi 
instead of \Qi\, etc.): 



where 



M 



jM _ 

In 



jM 

J- A — 



I 



2tt 2 



ln(l 



tM 



In 



2tt 2 D m D 4 D 5 P 2 QIQV"^ " Q<h " 
(M 2 + QI + QI- R M ) {Q\ - K?) {Ql - iff) 



1 + 



m 2 Q 2 Q\ 
r dft(Qi) dQ(Q 2 ) 



Sm 2 Q 2 Q 2 



7m ym 
/j PQ 1 L PQ 2 



J 2tt 2 2tt 2 D 4 D 5 


P 2 QiQ2 ' 


{Ql - iff) {Ql - iff) 




4m 4 Q?Ql 




r dn(Q ± ) dQ(Q 2 ) 1 


yM 7m 

^QxQi ^PQi 


J 2tt 2 2tt 2 D m D a 


PQ2QI 


{M 2 + QI + QI- R M ) 


{Ql - R?) 


Am 2 Q\Q 2 2 


1 


f dn{Q x ) dQ{Q 2 ) 1 


yM ym 
_ Z QiQ 2 Z PQ 2 



2tt 2 2tt 2 D m D 5 PQxQl 
(M 2 + g 2 + Ql - R M ) {Ql - iff) 

4m 2 g 2 g| 



7 2tt 2 2tt 2 






g?-i?r 






2m 2 g 2 ' 








1 


ym 
/j PQ2 


J 2tt 2 2tt 2 


^ 5 


PQ2' 


n2 nm 






2m 2 g| ' 






/■ dfi(Qx) rf^(g 2 ) 


1 


yM 
_ Zj QxQ2 



J 2tt 2 2tt 2 D m 
M 2 + Q 2 + Ql - R M 



Q1Q2 



R M = \l {m 2 + g 2 + g 2 ) 2 - 4g 2 g 2 , 
#™ = ^' + 4m 2 g 2 , i = 1,2. 



(5.9) 



(5.10) 



We shall also need the angular integrals where the massive propagator Dm is replaced by 
the massless one D 3 . We have set P 2 = —m 2 in the explicit expressions involving the square 
roots. Since the external momentum P flows only through the massive fermion propagators, 
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we do not need to deform the integration contour for the radial integrals over Q\ and Q\ (see 
the discussion in Refs. |3"5|). Note that I5 and Iq depend only on one variable, Q\ and 
Ql, respectively, whereas the expression for I 2 factorizes into a product of two single- variable 
functions. 

Other angular integrals that occur during the calculation, involving the factors D 3 / (D^D 5 ), 
Da/ (DmD§), D5/ {DmDa)i and D3/.D4, can be reduced to the basic ones, by noting that D 3 = 
QI + QI + 2Qi • Q 2 , L> 4 = Ql + 2P • Q u and D 5 = Q\ - IP ■ Q 2 , and with the help of the 
identities 



M 2 + Ql + Ql M 

Q1Q2 QlC 
m 2 + P 2 + Ql „ m 



PQ l 



( 7 m y 

In this way one obtains, for P 2 = —m 2 , 

dfi(Oi) dQ(Q 2 ) Q x -Q 2 1 (^p Ql ) 2 ( z pq 2 ) 



J PQ t 



1 . 



(5.11) 



2vr 2 2tt 2 




dn(Qi) dn(Q 2 ] 


p-Oi 


2w 2 2% 2 




dn(Q0 dn(Q 2 ] 




2tt 2 2tt 2 




dn(Qi) dn(Q 2 ) 


0i-0 2 


2tt 2 2tt 2 


P 4 



4 P 2 
QlQ 2 2 h-Qlh-Q 2 2 h + l 
Am 2 

1 { Z q[q 2 ) { Z PQ 2 ) 
4 Q2 

_[m 2 + g 2 + g 2 ]g 2 if + (m 2 + g 2 + g 2 ) if + Q 2 h-i 

AQl 

1 { Z QiQa) { Z PQi) 

4 o? 

[M 2 + Q 2 + Q 2 }Ql If - (M 2 + Ql + Ql) If -Qfh + 1 



40? 



(5.12) 



6 Two-dimensional integral representation for a^ byL;7r ° 

After having performed the angular integrations, the pion-exchange contribution to g — 2 can 
be written in a two-dimensional integral representation as follows: 



LbyL;7r 

v 



LbyL;7r°(l) 



LbyL;7r°(l) , LbyL;^°(2) 



(6.1) 



/•OO /"CO 

/ dOi / ^g 2 

JO 



^(Oi.Oa) / (1) (Oi,o 



13 



£ w gi (M Vi ,Q 1 ,Q 2 ) g$ v XQl,Q 



a 



LbyL;ir°(2) 



/ dQ x dQ 2 w g2 (M,Q l ,Q 2 )g ( $(Ql,Ql), 

-'0 JO » / f_»,r 



(6.2) 
(6.3) 



M=Mtt,M v 



with [see Eqs. Q-Q] 

fey 

(Qi,oi) 



9m v .(-QT, 
3 M 2 



J 7 ( Q2, 0) 



9m v . (0) 



My, ^ - Ml 



J- n-0-,*-,* ( Q 1; Q 



9m Vi (0) 
(i) 



(6.4) 



Note (see Table ffl) that in the WZW model we obtain gjjV (Qf, Q|) = 0, gf) (Ql,Q%) = 0, 
-Q\) = 0. On the other hand, for the VMD model we have the simplifications 



since g 



WZW 1 



fU{QlQl) = 0, gZSQlQt) = -QZiQlQH since f VMD (-Qj) = 0. 

The universal [for the class of form factors that have a representation of the type (4.6)1 
weight functions in Eqs. ( |6.2| ) and ( |6.3| ) are given by 



w fl {Qi,Q 2 ) 



7T 



-A{2m'-Qi)(Q{-Qi 



2\2 



Ql + Ml Qm 2 Q 1 Q 2 

(Ql + Q*-B!>){Ql-R?)(Q>-HZ 



x In 



1 + 



+ 



8m 2 Q 2 Q 2 

- Am 2 QlQ 2 2 - —M\Q\ + ^-MQl + Ql 



3QtQ 2 2 - Q'iQt + Q 



2/n4 



2m 4 



1 



R? 

Ql 



{Ql-QIYR 1 



1 



Ql 



+ — 2 Qi{Ql-Am i )RT 



2m 4 



QiiQ 



2 - 2m 2 ) 2 i? 2 m 



+ 



w gi (M,Q 1 ,Q 2 ) 



2m 4 

7T 2 



(^-w),^-w 0? - 2m ^, fl r«r 

1 



(6.5) 



Ql + Ml 6m 2 Q 1 Q 2 



4(2m 2 - Q 2 2 )(Qi - Q\ 



x In 



1 + 



(Ql + Ql-R°)(Qj-K?)(Ql-R™ 
8m 2 Q 2 Q 2 
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w g2 {M,Q 1 ,Q 2 ) 



+4(2m 2 - g 2 )[M 4 + (Ql - Q\f + 2M 2 (Qj + Q 2 2 )] 

(M 2 + QI + QI- R M )(Q 2 - R?)(Q 2 - if 



x In 



1 + 



8m 2 Q 2 Q! 

M 6 + 3M 4 Ql + 3M 2 Q\ + 3M 4 Q 2 2 + 2M 2 Q 2 Q 2 + 3M 2 Q 



M 2 



QiQ 



-(Ql-Q 2 2 ) 2 R° 



iff 

i - £) - ^OJW 



2m 2 )iff(l - |r 



+ 



M 4 + 2M 2 Q 2 + 2M 2 Q 2 + g 4 + g 4 - 2Q 2 Q 2 



6m 2 M 2 g!g 2 



|4[m 2 (g 2 -g 2 ) + 2g 2 g 2 ](g 2 -g 2 ) 

(g 2 + g 2 -i?°)(g 2 -^)(Ql-^) 
8m 2 g 2 g 2 

+4[M 4 m 2 + (Ql - Ql) \-1Q\Q\ + m 2 {Q\ - Q 2 2 )] 



x In 



x In 



1 + 



+2M 2 [g 2 g 2 + m 2 (g 2 + g 2 )] 

(M 2 + Q 2 + Q 2 - R M ){Q 2 - R?){Q 2 - iff) 



8m 2 g 2 g 

+M 4 g 2 + 2M 2 g 4 + M 4 g 2 - M 2 g 2 g 2 + 2M 2 g 4 

+ Ql(Qi + Ql)R°(l 

\ Vi/ 

iff 



pm 

+ Q^(Q^-3g?)i?°(i-^- 



Ql 



-Ql(M 2 + Q 2 l + Q 2 2 )R M (l - 



M 2 {M 2 + 2Ql + Q 2 2 )R™ 



(6.6) 



- g 2 (M 2 -3g 2 + g 2 ),R 

- M 2 iff iff 1 , 



pm 

^L) -M 2 (M 2 



3Ql + 2g 2 )iff 



(6.7) 



where 



R° = R 



M=0 



{Ql + QIY-AQIQI 



(6. 



Although the analytical expressions for the weight functions look quite complicated and 
involve terms with different signs and of different sizes, the sum of all terms leads to rather 
smooth functions of the two variables Q\ and Q 2 , as can be seen from the plots in Fig. ^| We 
have not shown the corresponding plots for w gi (Mv 2 , Qi, Q 2 ) and w g2 (Mv 2 , Qi, Q 2 ), since they 
look qualitatively similar. The contribution of the former is suppressed in the integral by the 
factor 1/My 2 in g$ v whereas the latter weight function scales as 1/M V2 for large My 2 . 
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Figure 3: The weight functions of Eqs. (|6.5|) - (|6~77|) . Note the different ranges of Qi in the 
subplots. The functions and w gi are positive definite and peaked in the region Qi ~ Q 2 ~ 
0.5 GeV. Note, however, the tail in Wf 1 in the Qi direction for Q 2 ~ 0.2 GeV. The functions 
Wg 2 (M 7T ,Qi,Q 2 ) and w g2 (Mv, Qi, Q2) take both signs, but their magnitudes remain small as 
compared to Wf^Qx, Q 2 ) and w gi (Mv, Qi, Q2)- We have used My = M p = 769 MeV. 
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The functions and w gi are positive and concentrated around momenta of the order of 
0.5 GeV. This behavior was already observed in Ref. by varying the upper bound of the 
integrals (an analogous analysis is contained in Ref. pOj ). Note, however, the tail in Wf 1 in the 
Qi direction for Q 2 ~ 0.2 GeV. On the other hand, the function w 92 has positive and negative 
contributions in that region, which will lead to a strong cancellation in the corresponding 
integrals, provided they are multiplied by a positive function composed of the form factors (see 
the numerical results below). 

As can be seen from the plots, and checked analytically, the weight functions vanish for 
small momenta. Therefore, as already noted before, the integrals are infrared finite. For large 
momenta the weight function Wf^Qi, Q 2 ) and w gi (M, Qi, Q2) have the following behavior: 

r (n n \ 2 [Qf - Q 2 (Ql - 2m 2 )i^] , „ ( 1 \ 
hm w fl (Q 1 ,Q 2 ) = 7T „, VO —3 , Q 2 fixed , 



lun^AQuQi) = - ^ + Qlfixed - 



2, 



13m 2 .„ / 1 



lim^W.Q) = ^_ + ^_j, (6.9) 
and 

2 M 2 [-18m 2 Ql + 5Qf + Qi(28m 2 - 5Q\)RT} 



\^ Wgi {M,Q u Q 2 ) = vr 9m2g| 

' 1 ' 

M 



+ 0[- f ), Q x fixed, 



22M 2 m 2 ( 1 \ 

= *^- + <%]- (610) 

Finally, the weight function w g2 (M,Qi,Q 2 ) behaves as follows: 
r (mrn n\ 2 [-Ql(Q 2 + m 2 ) + Q 2 (Q 2 - m 2 )R™} , n fl\ n , , 



-QUQi + Sm ^+QxiQj + Zm 2 )^] M 
3m 2 Q| + 1q| 



lim w g2 (M,Q 1 ,Q 2 ) = tt 2 L ^ IV ^ 1 ; 2 7,3 ^ + 75 , Qi fixed. 

lfim 4 / 1 \ 

hm^ 2 (M,Q,Q) = * 2 ^ + °[qe)- (6- 11 ) 

Since {Q11Q2) = in the WZW model, the corresponding integral a^ 1 ^ ^ involves 

only w f 1 and diverges as C In 2 A for some ultraviolet cutoff A . Varying the cutoff A in 
our numerical integration between 2 and 50 GeV, we obtain C ~ 0.025. Since all form factors 
tend to the WZW model for My — > 00, the results for a Lb y L;7r°(i) - n ^ e different models should 
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all scale as C In 2 My for large resonance masses, with the same coefficient C in front of the 
log-squared term as in the WZW model. We have numerically confirmed this observation for 
the VMD and LMD form factors. In fact, for the WZW form factor, the region Qi > Q2 in the 
integral produces the In 2 A term. Using the asymptotic expansion for large Q 1 of Wf 1 (Qi,Q 2 ) 
in Eq. ( |6.9| ), one can extract the coefficient of the leading-logarithm term from the integral, 
with the result 

( N c \ 2 (m 



c = 3 {wJ {tJ ■ < 612 » 

The corresponding value C = 0.0248 (for Nc = 3 and the input values for m and given in 
the next section) is in good agreement with the numerical determination quoted above. On the 
other hand, a ! f f yL ' ,7r0 must vanish as the lepton mass m tends to zero. Numerically, we find that 
it is indeed the case; see the next section. 

Before we present our estimate for the pion-exchange contribution to the muon anomalous 
magnetic moment in the next section, we would like to stress again that the two-dimensional 
integral representation ( |6.1| )-( |6~3| ) is valid for any form factor that can be written as shown in 
Eq. ( f4.6|) . If in the future experimental data for the off-shell form factor jF 7r o 7 * 7 * (— Q\, — Q 2 ) be- 
come available, the region below 1 GeV should also be measured with high precision. Provided 
the data can be fitted with a representation of the form factor belonging to the class we have 
discussed, Eqs. QS.lp-flPD can then be used to improve the present estimates. In this sense, 
our integral representation is similar to the familiar one that connects the vacuum polarization 
contribution to g — 2 and the experimental cross section a(e + e~ — > hadrons) |10[ . 



7 Numerical results 

The weight functions in Eqs. (|6.5|)-(|5~77|) are composed of rational functions, square roots 



and logarithms. The combinations of form factors from Eq. Q6.4J) that enter in the inte- 
grals are rational functions as well. It might therefore be possible to perform one or even 
both integrations over Q\ and Q2 analytically for some parts of our expressions.^ On the 
other hand, the numerical evaluation of the two-dimensional integral representation for a^^ 
in Eqs. ( |6.2| ) and ( |6.3| ) poses no real problems. The corresponding results for the differ- 
ent form factors discussed earlier can be found in Table |2|. Apart from the masses of the 
muon m = m M = 105.66 MeV and the pion M n o = 134.98 MeV, we have used the fol- 
lowing input values: a = 1/137.03599976, F n = 92.4 MeV, M Vl = M p = 769 MeV, and 
My 2 = M p i = 1465 MeV ||43|| . The errors from the numerical integration are much smaller 
than the last digits given in the table. 

As already announced in the Introduction we obtain a different sign for a^ byL;7r ° as compared 



to the latest two calculations [20, 21, 22 1. This is immediately visible from the plots of the 



weight functions and w gi which are positive and which are multiplied by positive functions 



/« and 

9m v i a ^ least for the WZW, VMD, and LMD form factors. We shall further discuss 
this important point in the next section. On the other hand the contribution aj^ byL;7r ^ is highly 

3 For the analytical integration it might be convenient to extract the factor Q1Q2 from the weight functions 
and to rewrite the integrals in terms of Q\ and Q\ . 
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Table 2: Results for the terms a^ byL;7r °W, a LbyL ;7 r°(2) anc j £ or pion-exchange contribution to 
the anomalous magnetic moment a^ byL;7r ° according to Eq. ( |6.1| ) for the different form factors 
considered. In the WZW model we used a cutoff of 1 GeV in the first contribution, whereas 
the second term is ultraviolet finite. In the LMD+V Ansatz we used hi = GeV 2 and h 5 = 
6.93 GeV 4 . 



Form factor 


a LbyL;^°(l) 


LbyL;7T°(2) 
a fJ, 


a iW> x 1Q io 


WZW 




0.095 


0.0020 


12.2 


VMD 




0.044 


0.0013 


5.6 


LMD 




0.057 


0.0014 


7.3 


LMD+V (h 2 


= -10 GeV 2 ) 


0.049 


0.0013 


6.3 


LMD+V (h 2 


= GeV 2 ) 


0.045 


0.0013 


5.8 


LMD+V \h 2 


= 10 GeV 2 ) 


0.041 


0.0013 


5.3 



suppressed and very stable with respect to the various form factors (with the obvious exception 
of the WZW case). As noted earlier, the corresponding integral converges even for a constant 
form factor as in the WZW model. As far as the numerics is concerned, the absolute value 



for ajf yL;,r for the VMD from factor is identical to the results quoted in Refs. [19], |2(| |2T], [Z| , 
whereas the one for the LMD form factor agrees, up to the sign, with the result for the form 
factor number 4 given in Ref . jjjj . 

Since the results for the cases of the LMD+V and VMD form factors are very similar, 
one might conclude that imposing the asymptotic l/Q 2 behavior of ^ r 7r o 7 * 7 *(— Q 2 , 0) for large 
spacelike momenta is more important than reproducing the QCD short-distance behavior for 
Q\ — > oo and Q\ — > oo, as does the LMD form factor. However, the slightly higher result 
for a^ byL;7r ° in the case of the LMD form factor might also be due to the fact that in this 
instance the slope of the form factor at the origin, as obtained, for instance, by the CELLO 
Collaboration f50fl , is not well reproduced, in contrast to the LMD+V and VMD cases (see the 
discussion in Ref. p5|| ). As noted in connection with the plots of the weight functions, it is 
mainly the region Q%% \ GeV that contributes in the integrals. 

For illustration, we decompose, in the cases of the LMD and LMD+V form factors, the 
total numerical result of Table according to the integrals corresponding to the different weight 
functions appearing in the sums in Eqs. (|6.2|) and (|6.3|), 



Lby L +>(i) = 0.015 + 0.042 + = 0.057 



LbyL;.0(2) = q.0016 — 0.0002 + = 0.0014 







a Lby L; ^(i) = 0.0026 + 0.0448- 0.0026 = 0.045 
a Lb yL; ^(2) = 0.0015 - 0.0002 - 1 x 10~ 6 = 0.0013 



LMD, (7.1) 
LMD+V (h 2 = 0) . (7.2) 



One observes that in both cases the term involving the weight function w 9l (My 1 , Qi, Q 2 ) gives 
the main contribution to the final result. 

For a fixed value of h 2 in the LMD+V form factor, our results are rather stable under 
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the variation of the other parameters. For instance, varying My i; My 2 , and h§ by ±20 MeV, 
±25 MeV, and ±0.5 GeV 4 , respectively, the result for aJf yL;7r ° changes by ±0.2 x 10~ 10 . On 
the other hand, if all other parameters are kept fixed, our result depends almost linearly on 
h 2 , at least for \h 2 \ < 20 GeV 2 . In this range aJf yL;7r ° changes by ±0.9 x 10~ 10 from the central 
value for h 2 = 0. As noted earlier, the experimental data for tt° — > e + e~ |42| , |43|j seem to favor 



slightly positive values for h 2 , although the bounds are rather loose. A better experimental 
determination of this decay rate would therefore be highly welcome. 

Thus, using the LMD±V form factor, our estimate for the muon anomalous magnetic mo- 
ment from the pion-pole contribution in hadronic light-by-light scattering reads 

where the error includes the variation of the parameters and the intrinsic model dependence. 

As far as the contribution to from the exchange of the rj or of the rf is concerned, a 
detailed short- distance analysis of the corresponding form factor in large- Ac QCD for nonzero 
quark masses as done for the pion is beyond the scope of the present work. Furthermore, as 
noted above, the integrals in Eqs. ( |6.2| ) and ( |673| ) do not seem to be very sensitive to the correct 
asymptotic behavior for large momenta. It is more important to have a good description at small 
and intermediate energies, e.g., by reproducing the slope of the form factor JF PS7 * 7 * (— Q 2 , 0), 
PS = rj, rj', at the origin (cf. the similar results obtained for the VMD and the LMD±V form 
factors). The CLEO Collaboration ]37[ has made a fit of the form factors jF r?7 » 7 » (-Q 2 , 0) 
and jF r? / 7 » 7 * (— Q 2 , 0), normalized to the corresponding experimental width T(PS — > 77), using a 
VMD Ansatz with an adjustable vector meson mass A PS . Taking their values = 774±29 MeV 
or Arf = 859 ± 28 MeV as the vector meson mass My in the expression of the VMD form factor, 
we obtain from Eq. ( |6.1|) 



< yL; "|vMD = +1.3 (0.1) x KT 10 , 

a L ; yUv '\ V MD = +1.2 (0.1) x 10- 10 . (7.4) 
The error reflects only the corresponding variation in A PS . Again, apart from the sign, these 



values are comparable with the numbers quoted in Refs. [21], |22], p9|1 . We do not expect the 
results with the LMD±V form factor to differ much from the numbers quoted above. Thus, 
adding up all contributions from pseudoscalar exchange, we obtain 

a L byL:P s = a L byL; .o + a ^v\ VMD + a ^'\ VMD = ±8.3 (1.2) x 10- 10 . (7.5) 

Finally, for illustration and as a check, the result for the pion exchange contribution to the 
anomalous magnetic moment of the electron reads 



LbyL;7T°(l) 



e 



2.0 x 10" 6 , 



a L byL ;.0(2) = 2 Q x 1Q -6 ; 

a iW = 5 j x 1Q -i4 ) (7 6) 

for the LMD±V form factor (with h 2 = 0). Despite the various factors 1/m 2 in the weight 
functions, the final result is much smaller than for the muon, as it should be. The results for the 
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other form factors are of similar size (LMD: a^ hyL ^° = 7.7 x 10" 14 , VMD: a L e hyL ^° = 2.6 x 1(T 14 ), 
except for the WZW case. Note that for the electron the second contribution a ^ h y L ^°( 2 ) [ s [ n 
general not smaller than the first one, in contrast to the muon case. 



8 Discussion and conclusions 

In this article, we have evaluated the correction to induced by the pion-pole contribution to 
hadronic light-by-light scattering. From the methodological point of view, the present approach 
differs from previous ones in two important points. First, we use a representation of the form 
factor .7 r 7r o 7 « 7 » which incorporates short- distance properties of QCD within an approximation 
to the large-A^ limit involving a finite number of resonances. It is important to stress that 
this approximation can in principle be improved by considering additional vector meson states, 
if additional information, subleading short-distance corrections, phenomenological and/or ex- 
perimental input, is provided. Second, for the corresponding class of form factors, we have 
performed the angular integration in the two-loop expression in an analytical way, using the 
method of Gegenbauer polynomials. This leads to an expression for aJf yL;,r in terms of weighted 
two-dimensional integrals over the moduli of the Euclidean loop momenta. These integrals were 
then evaluated numerically, and the results displayed in Table 2. 

Whereas the size of a^ byL;7r ° obtained in this way is quite comparable to the existing deter- 
minations, the sign comes out opposite. In view of the consequences implied by this result, it 
is clear that this discrepancy needs to be discussed. We first notice that our result is not an 
artifact due to the representation of the form factor we use. Indeed, repeating the same anal- 
ysis with the VMD form factor, we exactly reproduce the result quoted by previous authors, 
but again with the opposite global sign. Now, one might observe that the VMD form factor 
contributes only to the second term in Eq. ( |6.2|) , involving the positive definite (see Fig. ||) 
weight function w gi (Mv, Qi, Q2), so that our result in this case could be ascribed to the fact 
that we obtain the wrong sign for that function. In addition to the various checks that we have 
performed, we shall however put forward a few arguments that support the correctness of our 
result. In the case of the constant form factor F^SY*, only the first term in Eq. ( |6.2j ), involving 
the positive definite weight function w/^Qi, Q2), contributes. This integral is divergent, and 



behaves as |49], ^TJ 

/•A /-A / AT„ \ 2 

lim / dQ 1 / dQ 2 — WfAQiM = Cln ' A + (8- 1 ) 

A^ooJo JO \127T Z J H 7T / 

for a large ultraviolet cutoff A (the ellipsis stands for subleading In A divergences, and for finite 
terms). The coefficient C of this log-squared dependence can be computed analytically in a way 
[ KTJ that is largely independent of the methods used here. It is positive and agrees with the 
value we obtain from a numerical analysis of our formulas and from the asymptotic behavior for 
large momenta of the weight function Wf 1 (Qi,Q2)] see Eq. ( |b.l2| ). Next, as My becomes very 
large, the VMD form factor approaches the WZW form factor and therefore should reproduce 
the previous result, 

lim / dQ x dQ 2 w 9l (M v ,Q l ,Q 2 ) g^iQlQl) = C\n 2 M v + ••• (8.2) 
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with the same constant C as in Eq. (|8.1|) . We have checked numerically that this indeed happens 
with Eq. ( |6.2D and the expression of w 9l (Mv, Qi, Q2) given in Eq. (|6.4|). Thus, the signs of the 
contributions involving Wf 1 and w 9l must be correct. 

Finally, we mention that with the convention for the external momenta used in the third 
article quoted under Ref. [^], the right-hand side of Eq. (2.7) in that paper, which corresponds 



to our Eq. (|2.9|) , should come with a global minus sign. We also stress that our Eq. ( |2.9|) agrees 
with Eq. (2.9) of P2| . We have not been able to find a similar point of disagreement with Refs. 
|p0| f2~Tl . However, the lack of analytical expressions in intermediate steps precludes a more 
detailed comparison. Clearly, this matter needs to be investigated further in the future.^ 

Concluding this study, we find that the pseudoscalar-pole contribution to the hadronic 
light-by- light scattering amounts to 

a^ PS = +8.3 (1.2) x 10- 10 , (8.3) 

where we have used the LMD+V representation (fO|) of the form factor J r 7r o 7 . 7 « and the VMD 
form factor for the rj and r/, both supplemented with experimental information. The error 
given in Eq. (|8.3|) includes our estimate of the intrinsic model dependence and the variation 
of the model parameters. If we assume that the other contributions to the hadronic light-by- 
light scattering remain unchanged, our result implies that the difference between theory and 
experiment reduces to about 25 (16) x 10~ 10 , in the "least favorable" case, i.e., the vacuum po- 



larization analysis of Ref. [ II], which has the smallest value and the smallest error as compared 



to other recent evaluations ]T^, O, |I3|, [T5], [TEfl . We intend to use the same approach to study 
the full four-point function H^XpiQi, q2, #3) in order to obtain an evaluation of the complete 
contribution a L b y L ; had . We also plan to investigate in more detail the r\ and rj' pole contributions, 
although we do not expect the numerical values to change significantly as compared to the ones 



obtained in Eq. (7.4) with the VMD form factor. 
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Appendix 



This Appendix lists the properties of the form factor ^o^^^qf, q%) that are reproduced by the 
LMD and LMD+V representations. 

At the low-energy end, the form factor is normalized by the tt° — > 77 amplitude 

e 2 ^v> rr (0, 0) = A(n° -> 77) . (Al) 

In the chiral limit m q — > 0, q — u,d, s, this amplitude is fixed by the WZW anomaly to read 

For massive light quarks, this expression receives corrections. In particular, the pion decay 
constant in the chiral limit, F , is replaced by its physical counterpart F n — F [l + 0(m q )], 

- 77) = - T f^fr [1 + 0(m q )] . (A3) 



It turns out that the additional 0{m q ) corrections in this relation are numerically small |55 
so that one may drop them to a good approximation. 

The behavior of jF 7r o 7 * 7 * at short distances was recently studied in Ref. |2J| (see also the 
references therein) in the chiral limit. It turns out that the properties relevant for our present 
purposes still hold for massive light quarks, again upon replacing F by F n , 

Ljn W(AV, (V ~ M?) = § £ + ^ ^ + O (i) } . (A4) 

This expression holds up to 0(a s ) corrections. Note, however, that the Wilson coefficients 
corresponding to the two first terms of the short- distance expansion shown here are free of 
anomalous dimensions. 

In the large-iVc limit of QCD, the singularities of jF 7r o 7 . 7 * are restricted to single poles in 
each channel, 

2 2 ■ ( q^ q\ ) 



(q! - M^(qi - Ml) 



where the sum runs over the infinite tower of zero- width J PC = 1 vector resonances of large- 
Nc QCD. Beyond the fact that they must reproduce the short-distance behavior ( |A4j ) and be 
free of singularities, the functions Cy(qf, g|) are not further restricted by the (known) properties 
of large- Nc QCD. The representations J-^J^* and J-^ofRt v are obtained by truncation of the 
infinite sum ( |X5| ) to one, respectively two, vector resonances per channel. As noted in the text, 
as My — > 00, the form factor T^P^ t reduces to a constant given by the Wess-Zumino-Witten 

-LMD+V 
: cast; ui j~ 

is decoupled 



term. In the case of T^£t¥^ v one has to proceed in two steps. First the heavier resonance V% 



v 2 



Hm T^ V {ql ql) = Tgfc {ql q\) , (A6) 



23 



with ^ h h 

lim — r- = — cy, lim — — = 1 , lim — ^- = 0, i = 1,2. (A7) 

Then one lets — > oo as before. 
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